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I. INTRODUCTION 



In chemistry and atomic physics, 1 /N (or 1/ D) perturbation theory is very powerful in calculating energy eigenstates 
for many complicated systems [1-4] . The method is especially valuable when the applicability of Hartree-Fock methods 
is limited. In this study the behaviour of expansions is investigated at a very high order of approximation. For 
the sake of simplicity only radial potentials with spherical symmetry will be considered. The A^-dimensional non- 
relativistic Schrodinger equation for such potentials may be written as 
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where N is the number of spatial dimensions [5]. A brief review of this equation is presented in appendix B. 

The concept underlying 1/iV expansion methods is the use of perturbation theory on the expandable parameter 
A standard problem in perturbation theory would be a solvable Schrodinger equation plus a small perturbation, 
expressed as 



H = Ho+ gV 



(2) 



where Hq is the solvable Hamiltonian and gV is the small perturbation. In this notation g is a real expansion 
coefficient. To optain an approximation to the eigenstates of the perturbed Schrodinger equation one inserts a series 
for the wave function in powers of g as 
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and a series in energy as 
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Finally, these series are solved for the unknown functions u*^"^^ and the energy coefficients E'^'^\ 

For the standard problem this procedure provides an efficient way to obtain an estimate of the eigenstates and the 
energy of the physical system. Widespread application of this method is, however, limited since only problems with 
an expandable coefficient g can be solved. The dimensional dependency of the Schrodinger equation provides us with 
a hidden expandable variable of the potential. Roughly speaking, enlarging the dimensions of a physical system turns 
the system "classical" , and in the limit of large dimensionality a particle in a potential is fixed at the minimum; in 
some sense, h —i- 0. A classical analog could be a damped pendulum. A differential equation describing this is e.g.: 



mrO + bO + mgO =^ 



(5) 



where mrO <C hO (assuming small angles). In the case of a very heavily damped pendulum it is safe to disregard the 
Omr term. This implies that the second-order differential equation reduces to the first-order differential equation 



be + mgO = 



(6) 



Solutions of a first-order differential equation are unable to oscillate. The pendulum thus seeks the fixed point of the 
potential. 
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An A'^-dimensional particle, where N is very large, behaves somewhat like a pendulum in a highly viscous fluid. 
Its "damping term" (the first order differential term) in the Schrodinger equation is large, and hence its de Broglie 
wavelength becomes small. To leading order in 1/iV, the particle may be addressed by quasi-classical concepts. 
Accordingly, the fixed point of the potential provides a starting point for a power series expansion in 1/N for the 
energy eigenvalue and the wave function. 

Most l/7V-series provide accurate results even when only the first few terms are summed; it is thus more interesting 
to examine the behaviour of the series as the order of approximation increases. An interesting aspect of the 1/N- 
method is that it is perturbatively exact in the sense that there are no explicit approximations in the method. Most 
potentials lead to apparently asymptotic series at large orders of approximation. The approach here is to calculate 
the series to a significant order (30 - 100 terms) for miscellaneous potentials and then explore the behaviour of 
the series. This is accomplished by employing a computer algebraic program. Round-ofi^ errors are also taken account 
of in the calculations. 



II. THEORY 

How to obtain the l/A'^-series 

In order to obtain 1 /A^-serics for a given potential there are several ways to proceed. One procedure is to expand 
the wave function, the energy, and the potential in Taylor series [6]. This approach makes it easy to make an efficient 
computer algorithm. Before the Taylor series expansions are made, it is preferable to redefine the equations slightly. 
First, the A''-dimensional Schrodinger equation is transformed into 
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by redefining ijj{r) = r 2 4>{r). Now, k is defined by fc = A^ + 2L The energy eigenvalue is denoted by E'. 
The position r is then redefined in terms of k, r = Vkp or yr = p where y is The Schrodinger equation is 
subsequently recast into to a differential equation in p, as 
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If now V{r) = kV{p), then 
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In the large- A" limit one has an effective potential like 
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Its minimum is the energy kE^ which matches the energy at the minimum po- 

To ease the further derivations it is preferable to rescale distance as a; = Vk{p — po) and define '4^{p) = e^^^^ and 



V,s{x) = Ip{x)-' + V{p) - E(-'^ 
Introducing -0 in eq. (9) leads to a differential equation for U(x). 
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where W{x) = kVcs{x) + 
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. Now the actual Taylor series expansion may begin. Extending U{x), E, W{x) 



in Taylor series in x and y one obtains. 
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Coefficients of the type _Dq arc defined to be zero. The scries for the potential is extraordinary. For W{x) it turns 
out that for any power of x, say k {k > 1), one has solely terms of y in powers k ~ 2, k and k + 2. 



The recursion formulas 

Using the Taylor scries technique, scries expansions for the energy and the wave function have to be found by 
recursion. The easiest way to do this is to find universal recursion formulas [6] for the required coefficients and then 
use a computer to recursively calculate each coefficient one by one to a preferred number of terms. For the ground 
state, the recursion relations are established by the requirement that 
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From these equations it is apparent that 
and 

Likewise, for the coefiicients D'^ it follows that 
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For the coefficients one obtains 
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And for the energy one obtains, after some algebra 
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All coefficients of the type C", where j is negative, or where j is zero or negative, are assumed to be zero in the 
summation. 

Obvious extensions of this technique to include excited states are possible. We can always write 



i>={xy-A{y))e''^^^ 

for the first excited state, and 

^={x^y^+xyC{y) + B{y)) e^(-) 
for the second excited state, etc. Here, A{y), B{y) and C{y) are series expansions in y for the nodes in ij). 
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Subsequently, a differential equation for U [x) is obtained for the first excited state. 



{xy - A) {U{x) + U{xf - 2W{x) + 2E^ + 2yU{x) = 



Expressing 

in a new Taylor series 
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the previously obtained universal recursion formulas for the ground state are again useful. The recursion formulas for 
T and S are arrived at using the relation 
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One obtains 
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Once again 
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Successively for the second excited state it is found that 

(a;2y2 + Cxy + B) (u{x) + U{xf - 2W{x) + 2E^ +2y^ + 2 {2xy^ + yC) U{x) = 
Note that there is a misprint in [6], in that 2x^y should read 2xy^. Once again 
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Employing the Taylor expansion technique the sign of the coefficient has to be decided. The two signs lead to 
distinct dissimilar solutions. The choice leading to finite energy eigenvalues is the negative solution, and this is the 
one we must choose. 



III. THE RESULTS 



Exploying a computer algebraic algorithm programmed in MapleV [*], all coefficients were calculated to a selected 
order of approximation by recursion. 
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The energy estimates obtained using the method were compared with strict analytical or accurate numerical 
results from different sources [6-10]. In Table I, the results of the series for the harmonic oscillator and the Coulomb 
potential are displayed. The approximation to the appropriate energy eigenvalues is extremely good, indicating that 
in this case the Taylor expansion procedure is exceedingly reliable. The 1/A'' outcome for two additional potentials 
is also shown in the Table I. The results are excellent, but the energy sum starts to diverge at some point in 
the approximation (around 29th order or so), suggesting that the approximation is correct to about 6-7 digits. To 
illustrate the behaviour graphically, some plots of the partial sums of the energy approximation versus the order of 
approximation have been made for these potentials (see figures 1, 2 and 3). 

Some potentials, also dealt with in [6, 9], are recomputed here to a much higher order of approximation (see Tables 
II, III and IV). The results are very reasonable, and the high order of approximation make the results nearly exact. 
In particular, it is observed that the high-Z states approximate the correct energy quite well. In Table IV for the 
2nd excited states, only series which did not begin to diverge to around 30th order of approximation are presented. 
Most of the approximated partial energy sums for the 1st and 2nd excited states diverged at low I in an oscillatory 
manner around the correct energy eigenvalue. A high-Z quantum number entirely cancelled the divergence, leading 
to a seemingly convergent series. As seen from the Tables III and IV the results for the harmonic oscillator and the 
Coulomb potential are excellent and correspond exactly to the analytically known energy eigenvalues. 

The surveyed energy series are frequently seen to oscillate around a certain energy value when plotted against the 
order of approximation. The exact energy eigenvalue is found to be around this energy. The series begins to diverge 
at some order of approximation but still varies around the same definite energy. The oscillating series are often seen 
to fluctuate not termwise, but with two terms or more above the specific energy value and then two terms or more 
below, etc. 

Some convention is required in order to discuss the results when this behaviour is seen. In order to discover the best 
energy eigenvalue, the last (first) term above and first (last) term below (or opposite) the center of oscillation with 
least variation across it is found. The result presented in the Tables is then these two terms. For a normal oscillating 
series, the two terms with least variation across the center of oscillation are presented. In Tables II, III and IV the 
normal oscillating behaviour is observed. The order of approximation of the partial sums employed is noted in the 
Tables. For the convergent series, only the highest order term is presented. 

To examine the oscillations more carefully another potential also discussed in [6, 9], namely the linear potential 
22r, will be analysed, (the curious value of the constant is chosen such that the results in this paper can be directly 
compared with those in [6, 9]) The results are presented in Table V. For this potential the energy eigenvalue is 
calculated to about the 29th order using the 1 /N method for the first three excited states. The results show partial 
energy sums for the ground, the 1st and 2nd excited state seemingly swinging around the correct energy eigenvalue. 
Some graphs showing the oscillations of the partial energy sums are presented (see figures 4, 5 and 6). 

Different radial potentials for which "exact" results from [9] arc used, arc shown in Tables VI. For most of these, 
satisfactory agreement with exact values is obtained. Many of them are not termwise oscillating series. The order of 
approximation employed gives a hint about the convergence of the expansion. The r'' potentials where A; > 3 are seen 
to diverge very quickly. An empirical formulation of this principle would be "the higher the power of the potential, 
the poorer the convergence of the resulting 1/A^ energy series." 

In order to investigate the applicability of the 1/A^ procedure for potentials with analytically precise eigenvalues, the 
method is now applied in the case of fabricated potentials where the potential is set up from the Schrodinger equation 
(see appendix A). These examples illustrate the reliability of 1/A^ progressions in the case of explicit, analytically 
known solutions and are hence interesting for our purposes. The results in Table VII are obtained using potentials 
calculated from a selected eigenfunction. Most of the constructed potentials are seen in particular to have a not 
termwise oscillatory behaviour. 

Miscellaneous "quark" potentials and a particular double- well potential have also been developed in 1/A^ scries (see 
Tables VIII and IX). The double- well potential has also been handled in [10] using ordinary perturbation calculations. 

These results are not quite as good as those encountered for the other potentials, but the Taylor method still yields 
a fairly good approximation to the proper eigenvalue. 

The double-well potential ( ^ ) where i? is a constant) is in fact a rather peculiar potential because in a normal 

perturbation calculation one has to add powers of small terms proportional to ~ exp(— |i?^) in order to take into 
account quantum-mechanical tunnelling [10] in the perturbation series. Such tunnelling considerations have not been 
included here, but clearly for certain potentials such terms proportional to e.g. exp(— aiV) should have an effect on 
the approximated energy eigenvalues. Maybe that is why the result for the double- well potential deviates slightly 
from the correct energy eigenvalue. This should of course be investigated in more details before using the Taylor 
expansion method more generally. To demonstrate the behaviour of some of the approximations, plots of the partial 
energy sums versus the order of convergence for some of the potentials are presented (see figures 7, 8, 9, 10, 11 and 
12). 

Most series are generated to about the 29th order of approximation with a numerical precision of about 100 digits. 
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In fact, it is feasible to generate the series to the 100th (or higher) order of approximation with a numerical precision 
of about 1000 digits. Further, it is conceivable to calculate the l/7V-series exactly using calculations with fractions 
and fractional roots to about the IGOth order of approximation. The most important potentials are the harmonic 
oscillator, the Coulomb potential, and with the constructed potentials examples where analytic exact solutions are 
available. 

Shanks resummed sums could have been used to resum the series. In a Shanks resummation Shanks extrapolants 
are calculated from the partial sums, as 

p p p2 

Pn+1 + Pn-1 ^ '^Pn 

where Sn is the nth Shanks extrapolant and P„ the nth partial sum. See [11]. In this paper it has been chosen not 
to do so because this question is an auxiliary aspect of the technique. 



IV. DISCUSSION 



The major aspect of the long term behaviour of the l/iV-series is of course their apparently asymptotic behaviour. 
As demostrated, many 1/A'' series match the exact energy satisfactorily when only a certain number of terms are 
summed, but as further terms are added the sum diverges oscillatingly to infinity. The plots presented clearly show 
this behaviour with some minor differences. Convergent 1/N series are also seen, such as the Coulomb and the 
harmonic oscillator potentials series, together with the high-Z states 1/N series (i.e. 1/k). This behaviour is very 
similar to that of asymptotic series. Most perturbation series are actually asymptotic. An asymptotic series is, briefly 
stated, a series which begins to converge towards a finite value, but which in the long run diverges. Mathematically, 
an asymptotic power series in - (here x plays the rule of N) is a series for which 
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for all zero and positive n, (Poincare's definition) [12] or [13] where f{x) is a function and the sum is a partial sum 
for the asymptotic power series of f{x). 

If the 1/A^ series arc asymptotic, the Coulomb and harmonic oscillator potentials woifld then just be special cases 
in which the divergence is extremely slow. Comparing the eigenfunctions for the Coulomb and harmonic oscillator 
potentials, they are both of the type exp(r'^), with A: = 1 and = 2, respectively. As a numerical experiment in this 
paper, eigenfunctions of the type exp(r'^) are used to find the corresponding potential with eigenvalue fixed ai E = \ 
(see appendix A). These potentials arc then expanded using the 1/A^ method, (see Table VII). Examining this table 
one notices that the nearer k is to 1 or 2, the later the divergence of the resulting series, suggesting that the point of 
divergence is related to k. For instance, with k equal to 1.2 or 1.15, the convergence is good to 25-30 orders or so. 

One possible hint could be the radius of convergence for the Taylor series of the U{x) function, but further inves- 
tigations examining the 1/A'^ series for, e.g., U{x) = cosh(.T) have shown that there is no simple relationship between 
the radius of convergence of the Taylor series of the U {x) function and the point of divergence of its corresponding 
1/A'' series. Furthermore, for A; = 3 the resulting series is divergent, like the series for k between 1 and 2. In a more 
comprehensive treatment this should be inspected more carefully, and the apparent relation between k and the order 
of convergence of the resulting 1/A^ series surveyed. 

The series reported in this article are derived to a very high numerical accuracy. Strong evidence that numerical 
round-off errors have not affected the results have been presented. By calculating results for different power-law 
potentials it has been found that results obtained "exactly" (by fractions and fractional roots) and results obtained 
with an exactness of 1000 digits and 100 digits, respectively, are nearly identical, and that differences occur only at 
about the 60th- 70th digit. 

The purpose of this article has been to consider the convergence of quantum mechanical 1/N series in detail. Many 
of the energy series presented arc extremely close to the exact energy eigenvalues. The results clearly demonstrate that 
1/A'^ methods can be used to obtain approximate energy eigenvalues for different physical and chemical potentials. 

The method has proved to yield fine results not only for the ground state but also for the 1st and 2nd excited 
states. For the harmonic oscillator and Coulomb potentials the results are extremely accurate. 

The question of the convergence of the series has been settled, and much evidence points towards the 1/A^ series 
being asymptotic. Some of the series are so rapidly convergent at low orders that it will be very difficult to see that 
they diverge at all. 
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It has been established that numerical round-off errors are completely insignificant here. This is a very important 
point because it demonstrates that the divergence problem is a feature of the method itself and not a product of the 
computer arithmetic calculations. 
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APPENDIX A: HOW TO CONSTRUCT THE POTENTIAL FROM A GIVEN GROUND-STATE 

WAVEFUNCTION 



If (f){r) = exp(^(r)) is an eigenfunction of 

l_(f_ _ l N-1 d 
2rfr2 ~ 2 r dr 

with eigenvalue E, then the potential must be 
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APPENDIX B: THE SCHRODINGER EQUATION IN N DIMENSIONS 

In this appendix we shall see how to derive the most important equation in this paper, namely the radial Schrodinger 

equation in dimensions (sec also [5]). 

The Laplace operator in polar coordinates can be written [5] as 
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where eo = r,h = Uf^^ hj and /if = X^^i (^) • 

The expressions for the spatial coordinates in N dimensions are given by 
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X2 = r sin 9i sin 62 sin 63 • • • sin 1 



Xi = r cos 6i-i sin 6i sin 9i+i ■■ -sin On -1 



xn = rcos^jv-i 
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for N >3, {xi = r cos6»i, X2 = rsinOi; N = 2) and where < r < 00; < 6I1 < 27r; < < tt; 2 < i < iV - 1 and 
The functions hj are given by the following expressions 



ho = 1 

hi = r sin^2 sin^s • • • sin^jv-i 

hj = r sin 9j+i sin 9j+2 • • • sin 0]v_i 



Inserting the expression for hj in equation (Bl) leads to 
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(B4) 



10 



In general, the Hamiltonian can be written [m = 1 = H) 



1 «2 f 2 

2 2 2r2 



(B5) 



separating angular momentum and radial momentum. Hence, from equation (B5) and equation (B4) one can define 
total angular momentum operators as 
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In three dimensions, Lf is simply Li, and L| is L^. Also, it is seen that has the same expression independently of 
the number of spatial dimensions. 
Equation (B4) may now be recast as 
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The eigenvalues of the operator L'j^_-^ will subsequently be determined. 

Only total angular momentum operators have so far been addressed, but as in three dimensions, angular momentum 
components also exist. They are written as 
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where pj is defined as 
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Considering the potential combinations, there are distinct (not just differing by a sign) operators of the 

pattern Lij . The generalized angular momentum operators can then be obtained in polar coordinates from equations 
(B2), (B9) and (B3). The results, i.e. the operators Lij, are then as follows 
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L2j = +i 
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for iV > 3, {ho — l,hi = r;N = 2) [5]. For = 3, it is noted that L12 is simply Lz, -£13 is —Ly and that L23 is Lx- 
As equation (B6) shows , L|, . . . , clearly commute because they depend on distinct angles, meaning that they 

have simultaneous eigenfunctions. Since is a sum of squares of Hermitian operators of the form {xiPj — xjPi), it 
is known that the eigenvalues of will be real and non-negative. Furthermore, the eigenfunctions arc orthogonal 
because they are eigenfunctions of Hermitian operators. Assuming a proper normalisation, they may be written as 
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Y (Ajv-i, . . . , Ai) = n^_i,...,Ai (^1, . • • , ^iv-i) = n ®' ^'-1) 
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where Qi (A;, Ai_i) is a function only of 9i, and 61 (Ai, Aq) = 0i (Ai) and A, is the associated eigenvalue of L?. 
Clearly, Qi (Aj, Ai_i) will satisfy 

L?ei(Ai) = Aiei(Ai) 



Li (Ai_i) Qi {Xi, — XiQi (Aj, Ai_i) 



(j-'i i^i-i) ~ Lf^i (Aj_2)^ Q'i-l {Xi-i, Xi-2) Q'i (Xi, Ai_i) 

= (Aj — Ai_i) Gi_i (Ai_i, Ai_2) Qi {Xi, Ai_i) 



3 = 1 



for 2 < i < iV - 1. 
Since 



(B12) 



(B13) 



the eigenvalues of \^Lf (Ai_i) — Lf_-^ (Ai_2)j = Aj — Ai_i must be non-negative and real simply because the sum is 
the diagonal matrix elements of the Hermitian operator {xjPi+i — Xi+ipj)"^ [5]. Clearly 

Ai > Ai_i (B14) 

which implies that 

Aiv-i > Ajv-2 > ••• > Ai > (B15) 

because Ai = > 

At N = 3,Xi=ll= 777,2 x2=l2{h + l) = l{l+l)=l{l + '2-l) 
This suggests that in (z — 1) dimensions 



Ai = li{k + i~ 1) 



(B16) 



Although the value of k is unknown, a sensible guess would be k € No since it should be legitimate for N = 3. By 
the principle of induction it is feasible to prove this proposition by proving that if 



ih-2) Qi-i {h-i, li-2) = k-i {h-i +i-2) Gi-i (^i-i, ;i-2) 



(BIT) 
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where k-i G No, 
then 



(B18) 



/j e No, and for a given Zj = 0, 1, 2, .. one has Zi_i = 0, 1, 2, .., Zj. 
Proof: 

Using (B6), L?(ii_i) may be transformed into 



ih-i) 



If ladder operators are defined as 



(z- l)cot0i— r-T^ 

odi sm Oi 



d 



{h-i) = - k-icotOi 
d 

L. = -— -(Zi_i+i-2)cot6»i 



then Lf (/i-i) can be written as 



Ll = L+ - 1) LT {k_,) + {k_, +t^2) - 1) 

= L- {k^i + 1) L+ + Z,_i + ^ - 1) 

By multiplying (B21) from right and left with {k-i) one finds 

+ 1) Lf = L+ ^ 

- 1) LT = ir Li 

If now (B22) operates on Oj (Aj, U-i) it leads to 

L2(/._i + i)i+(z._i)e,(Ai,/i_i) = Ai4+(/i_i)ei(Ai,/i_i) 

= Ai^(Ai,Zi_i)ei(Ai,/i_i + 1) 
L?(/i_i-l)Lr(Z._i)ei(Ai,Zi_i) = A,L-(Z,_i)e,(A„Z,_i) 

= Aj-B (A,, 6j (Aj, Zj-i — 1) 



(B19) 



(B20) 



(B21) 



(B22) 



(B23) 



where A, B arc normalizing factors. It turns out that B (A,, + 1) = A* (Aj, U-i) [5]. 
Thus, from (B21) 

-t'i (^i-l) ©i (Aj, Zj-l) = 

= {\A iX,,k-i) f + k-i ik-i +1-1)) Qi {Xi, h-i) 
=^ Xi - k-i {k-i + i - 1) = \A{Xi,li-i)\^>0 (B24) 

As a consequence of (B15), i.e. Aj > 0, there must exist a maximum value for li-i, say li, for a given Aj, such that 

L+{k)@i{Xi,li) = (B25) 

where k € No, because li-i G No. This leads to 

Xi =k{k+i- 1) (B26) 

where li is a positive integer or zero. 

To see what the possible values of are, given li, one notes that clearly = li is plausible and that (B23) implies 
that 



L-{li)Q, iX,,k) 



(B27) 



is an cigenfunction of Lf (li — 1) with eigenvalue A^. Hence, U-i = — 1 is viable. Continuing this argument and 
introducing Zj_i = — 1 into (B23), one discovers that L? {k — 2) has the eigenvalue Aj. Consequently, another 
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feasible value for is Zj — 2. Since is the lowest probable value of U-i by presumption and one diminish Zj_i one 
by one, it is unquestionably seen that the possible values of li^i are li, U — l,li — 2, .., for li € Nq. 
Prom (B15) it is demonstrated that Xn-i, and accordingly In-i, have no upper boundaries. It thus follows that 
In-1 = 0, 1,2, . . ., and in addition it is seen that I2 {h + 1) > Zf > 0, leading to h < \h\- Since (B17) is believed to 
be correct for i = 3, it will also be valid for i = 4, and therefore by induction 



L2e,(z„z,_i) = kik + t-i)ei{k,k-i) 



(B28) 



where l<i<N — l,0<l2<l3<U<---< Zjv-3 < In-2 < In-i < 00 and li < IZ2I 
Equation (Bll) may now be written as 



and 



N-l 



Y{lM-i,...,h) = n ^iih,h-i) 



i=l 



(B29) 



Qn-1, . . . . . . = + i — 1)Y (Zjv-1, ■ ■ ■ ■ ■ ■ ,h) 



(B30) 



Substituting i = N — 1 into (8), one finally arrives at 



"2rfr2 



l N-1 d l{l + N-2) 
2 r dr^ 2r2 



cj>{r) = Ecj>{r) 



(B31) 



as the A''-dimensional radial Schrodinger equation. 



[*] MaplcV, V. 3.0, Maple and MaplcV are trademarks of Waterloo Maple Software 

[1] M. Dunn, T. C. Gcrmann, D. Z. Goodson, C. A. Traynor, J. D. Morgan III, D. K. Watson, and D. R. Herschbach, J. 

Chem, Phys. 101, 5987 (1994). 
[2] D. K. Watson and D. Z. Goodson, Phys. Rev. A 51, R5 (1995). 

[3] Shi- Wei Huang, David Z. Goodson, Mario Lopez-Cabrera, Timothy C. Germann, Phys. Rev. A 58, 250 (1998) 

[4] V. S. Popov and A. V. Sergccv, Phys. Lett. A 193, 165 (1994). 

[5] James D. Louck, J. Mol. Spcctrocopy 4, 298-333 (1960). 

[6] Leonard D. Mlodinow, Michael P. Shatz. ,J. Math. Phys. 25(4), 943 (1984). 

[7] S. M. Blinder, ,J. Math. Phys. 25(4), 905 (1984). 

[8] Barnana Roy, Rajkumar Roychoudhury, Pinaki Roy, J. Phys. A21, 1579-1588 (1988). 
[9] T. Imbo, A. Pagnamenta, U. Sukhatme, Physical Review D, Vol. 29, No. 8 (1984). 
[10] R. Damburg, R. Propin, V. Marty shchenko, "Large-order perturbation theory for the 0(2) anharmonic oscillator with 
negative anharmonicity and for the double-well potential" in "Large-Order Behaviour of Perturbation Theory" , Edited by 
J. C. Lc Guillou, .L Zinn- Justin, Vol. 7, (1990) North-Holland, Elsevier Science Publishers. 
[11] Carl M. Bender, L. D. Mlodinow, N. Papanicolaou. Physical Review A, Vol. 25, No. 3 (1982). 

[12] R. B. Dingle, "Asymptotic Expansions: Their Derivation and Interpretation", Academic Press Inc, (1973) London and 
New York. 

[13] E. T. Copson, "55 asymptotic expansions", Cambridge University Press, (1967) printed in Great Britain. 



14 



-0.5010 
^ -0.5008 
5 -0.5006 
II -0.5004 
- -0.5002 

II -0.5000 ^ 

J -0.4998 T 
^ -0.4996 / 
5 -0.4994 ] 
H -0.4992 ■ 

-0.4990 ' -J 

10 20 30 

Order of approximation 



FIG. 1: Long-term behaviour of the 1/iV-energy series for the Coulomb potential. The dotted line represents the exact energy 
eigenvalue. 
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FIG. 2: Caption as for figure 1 for the harmonic oscillator potential. 
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FIG. 3: Long-term behaviour for an analytic solvable supersymmetic potential. Pot: + j^^^, t = 0.1 and g = 0.1. Angular 
momentum zero. [8] 
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FIG. 4: Graph of the energy series for the ground state of the linear potential: 25 r with zero angular momentum: / = 0. 

It is seen that the convergence for the first 10-20 terms is excellent. The series starts to diverge oscillatingly at about the 25th 

order of approximation. 
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FIG. 5: Caption as for figure 4 for the 1st excited state. 
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Order of approximation 



FIG. 7: Energy series for the ground state of the potential with zero angular momentum as a function of the order of 
approximation. As seen the convergence of the series is good up to a certain point, but at about the 20th order of approximation 
it starts to diverge. The series is seen to be oscillatory. 
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FIG. 8: Caption as for figure 7 for the potential r^. Now the divergence is much faster, and the series diverges even before it 
has reached its exact value (the power effect). 
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FIG. 9: Caption as for figure 7 for the In(r') potential. 
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FIG. 10: The ground state of the quark potential (^2.34)'^ ~ with zero angular momentum is seen to diverge at approximately 
the 20th order of approximation. 
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FIG. 11: Caption as for figure 10 for the quark potential +r*. The divergence just starts much faster than in figure 10. 
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FIG. 12: Energy series for the ground state of the double-well potential — g^^; ' , R = 4 with zero angular momentum. It is 
seen that the 1/N series converges to about the 20th order of approximation. It is also seen that oscillatory divergence does 
not occur for this potential. 
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TABLE I: Table showing the convergence of the 1/N series for 4 different potentials (2m = 1 = S), except for the Coulomb 
potential in which (m = 1 = K). The exact energy for the last potential is calculated in [8] with a supersymmetric method. 



The partial sums of the energy approximation 



^2 I 10000 ' 



1+1 



Order of approximation 



2.1213203433 
2.9526299690 
3.0036770689 
2.9997496105 
3.0000010566 
2.9999982706 
3.0000000000 

3exact \J\ 



-0.22222222222 
-0.49108367627 
-0.49993508216 
-0.49999961670 
-0.49999999794 
-0.50000000000 
-0.50000000000 



2.2154999 
3.0712441 
3.1238941 
3.1198959 
3.1200487 
3.1200752 
diverges 



Exact results for the eigenenergies 
-0.5exact[7] 3.120081 [8] 



2.94 



2.07828211 
2.89354432 
2.94360413 
2.93982608 
2.93996827 
2.93999831 
diverges 

■exact super sym 



1 

5 

10 
15 
20 
29 
100 



. [8] 



TABLE II: Results for the ground state energy (n = 0), AT = 3, for different potentials {m = h = 1). The "exact" results refer 
to [6, 9]; the calculated results are referred to as Sxayior- 



Potential I ^^shiftcd method [9] i^Tayior ^^exactnum [6, 9] Uscd Order 

_2i. 7^-0.2 Q -2.68601 -2.685882 -2.686 29 

_2i.7^-o.2 ^ -2.34494 -2.344946 -2.345 29 

_2i.7^-o.2 2 -2.15626 -2.156260 -2.156 29 

_2i. 7^-0.2 3 -2.02906 -2.029065 -2.029 29 

_20.8^-o.8 Q -1.21870 -1.218693 -1.218 29 

_20.8^-o.8 J -0.50044 -0.5004397 -0.500 29 

_20.8^-o.8 2 -0.29470 -0.2946959 -0.295 29 

_20.8^-o.8 3 -0.20191 -0.2019137 -0.202 29 



TABLE III: Results for the 1st excited state {n = 1), N = i, for different potentials {m = h = 1). The "exact" results refer to 
[6, 7, 9]; the calculated results are referred to as Sxayior- 



Potential 


/ 


-Bshifted method [9] 


£'Taylor 


£;exact„„„ [6, 9] 


Used order 


_2l. 7,^.-0.2 





-2.25483 


-2.253515 - (-2.25314) 


-2.253 


16-17 


_2l.7^-0.2 


1 


-2.10103 


-2.100738 


-2.101 


29 


_2l.7^-0.2 


2 


-1.99015 


-1.990056 


-1.990 


29 


_2l.7^-0.2 


3 


-1.90491 


-1.904867 


-1.905 


29 


_20.8^-0.8 





-0.46282 


-0.462291 


-0.462 
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1 


-0.28071 


-0.280648 


-0.281 


29 


_20.8^-0.8 


2 


-0.19493 


-0.194912 


-0.195 


29 


_20.8^-0.8 


3 


-0.14635 


-0.146342 


-0.146 


29 


Potential 


I 




-^Taylor 


-Boxact [7] 


Used order 


_ 1 

r 







-0.12500000 


-0.125exact 


29 


_ 1 

r 


1 




-0.05555556 


-0.05exact 


29 


_ 1 

r 


2 




-0.03125000 


-0.03125exact 


29 


_ 1 

r 


3 




-0.02000000 


-0.02exact 


29 


_ 1 

r 


4 




-0.01388889 


-0.0138exact 


29 


V 







3.50000000 


3.5exact 


29 




1 




4.50000000 


4.5exact 


29 


2' 


2 




5.50000000 


5.5exact 


29 


2' 


3 




6.50000000 


6.5exact 


29 


2' 


4 




7.50000000 


7.5exact 


29 
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TABLE IV: Results for the 2nd excited state (n = 2), N = 3, for different potentials (m = h = 1). The "exact" results refer 
to [6, 7, 9]; the calculated results are referred to as ^Taylor- 



Potential 


I 


-Eahifted method [9] 


-E-Taylor 


i;exact„„„ [6, 9] 


Used order 


_2l.7^-0.2 


1 


-1.95147 


-1.950722 


-1.951 


29 


_2l.7^-0.2 


2 


-1.87535 


-1.875032 


-1.875 


29 


_2l.7^-0.2 


3 


-1.81266 


-1.812502 


-1.812 


29 


_20.8^-0.8 


1 


-0.18745 


-0.187318 - (-0.187320) 


-0.187 


28-29 


_20.8j,-0.8 


2 


-0.14202 


-0.141980 


-0.142 


29 


_20.8^-0.8 


3 


-0.11280 


-0.112788 


-0.113 


29 


Potential 


/ 




-E'Taylor 


-Eexact [7] 


Used order 


_ 1 

r 







diverges 


-0.05exact 


29 


_ 1 

r 


1 




-0.03125000 


-0.03125exact 


29 


_ 1 

r 


2 




-0.02000000 


-0.02exact 


29 


_ 1 

r 


3 




-0.01388889 


-0.0138exact 


29 


_ 1 

r 


4 




-0.01020408 


-Mexact --0-01020408 

tJvj exact 


29 


2' 







5.50000000 


5.5exact 


29 


2' 


1 




6.50000000 


6.5exact 


29 


2' 


2 




7.50000000 


7.5exact 


29 


2' 


3 




8.50000000 


8.5exact 


29 


2' 


4 




9.50000000 


9.5exact 


29 



TABLE V: Table showing partial energy sums for the potential: 2 2r. Results for the first three excited states are shown 
{m = l = h). 



Ground-state energy Ist excited state energy 2ii(l excited state energy Order of approximation 

9.35240 16.35180 22.70413 10 

9.35229 16.35175 21.52097 15 

9.35153 16.35274 22.62047 20 

9.38745 16.44861 21.61422 25 

9.54245 16.35286 22.45952 26 

9.24170 15.68014 21.13008 27 

7.95750 15.92972 23.45280 28 

9.26333 21.51171 25.79375 29 

Exact energies [6, 9] 

9.35243exact 16.3518exact 22.08224exact ~ 

Shifted method [9] 

9.35243 16.32636 22.02319 ~ 
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TABLE VI: Results for the ground state energy (Z = n = 0), iV = 3, for different potentials (2m = h = 1). For the "exact" 
results, reference is made to [9]; ^Taylor refers to the calculated results. 



Potential 


^'shifted method [9] 


^-Taylor 


Sexactnum [9] 


Used order 




-0.29888 


-0.29880 - (-0.29931) 


-0.29609 


27-28 




-0.22035 


-0.22038 


-0.22029 


29 


0.15 

r 


1.32795 


1.32781 - 1.32797 


1.32795 


15-16 


^0.5 


1.83341 


1.83361 - 1.83287 


1.83339 


13-14 


^0.75 


2.10815 


2.10769 - 2.11223 


2.10814 


15-16 




2.70806 


2.70780 - 2.71299 


2.70809 


14-15 




3.45111 


3.46173 - 3.43840 


3.45056 


11-12 




3.80139 


3.81467 - 3.37768 


3.79967 


7-8 




4.09146 


4.19331 - 4.03767 


4.08916 


6-7 


Iu(/0 




1.01157 - l.Oilli 


1.0 il:! 





TABLE VII: Results for the ground state energy with (I = n = G) and iV = 3, for different constructed potentials with 
(2m = h — 1). The exact energies are chosen to be 1; Sxayior refers to the calculated results. 



Eigenfunction 


-^Taylor 


Used order 


exp(-r»-»") 


1.0005 - 0.99917 


12-13 


exp(-r-°-«^) 


1.000005 - 0.999993 


29-30 


exp(-r°-^°) 


1.00003 - 0.999998 


23-24 


exp(-r«-95) 


1.0000 - 0.99999 


24-25 


exp(-r^'"-') 


1.0000 - 1.00000 


24-25 


exp(-rii°) 


1.0000 - 0.99998 


25-26 


exp(-rii5) 


1.0001 - 1.00000 


18-19 


exp(-ri-20) 


1.0000 - 0.99998 


22-23 


exp(-ri-25) 


1.0000 - 0.99995 


18-19 


exp(-ri-3°) 


1.0000 - 0.99991 


18-19 


exp(-ri-35) 


1.0001 - 0.99985 


18-19 




1.0002 - 0.99975 


18-19 


exp(-ri-*5) 


1.0004 - 0.99960 


18-19 


exp(-ri-5°) 


1.0004 - 0.99953 


12-13 


exp(-ri-55) 


1.0005 - 0.99958 


12-13 


exp(-ri-«") 


1.0006 - 0.99960 


12-13 


cxp(-r^-''5) 


1.0008 - 0.99960 


12-13 


exp(-r^™) 


1.0002 - 0.99996 


15-16 


exp(-ri''«) 


1.0005 - 0.99990 


15-16 


cxp(-ri-»0) 


1.0010 - 0.99973 


15-16 


exp(-ri«5) 


1.0015 - 0.99944 


14-15 


exp(-ri-»°) 


1.0019 - 0.99963 


14-15 




i.OOKi - 0.<)<)<)(i8 


ii-io 



TABLE VIII: Results for the ground-state energy (/ = n = 0), A'' = 3, for different potentials (2m = h = 1). For the "exact" 
results, reference is made to [9], again Sxayior refers to the calculated results. 



Potential 


£'Bhifted method [9] 


i'Taylor 


£'exact„um [9] 


Used order 


6.8698r°-^ - 8.064 


-0.31914 


-0.31936 - (-0.31883) 


-0.31917 


11-10 


(2.34)^ 0.52^ 


0.47880 


0.48049 - 0.47990 


0.47811 


10-9 




4.65061 


4.66212 - 4.64267 


4.64881 


11-10 
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TABLE IX: Double- well potential: gj^i , -R = 4. Results for the ground-state energy [l = n = 0), {m = I = h), N = 3. 



i^num [10] 


-^Taylor 


-^perturbation [10] 


Used order 


0.483053433 


0.483018 - 0.483015 


0.483053390 


13-12 



